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We have extensively mapped a sample of dense molecular clouds (L1512, TMC-1C, L1262, Per 7, L1389, 
L1251E) in lines of HC 3 N, CH 3 OH, SO and C 18 0. We demonstrate that a high degree of chemical 
differentiation is present in all of the observed clouds. We analyse the molecular maps for each cloud, 
demonstrating a systematic chemical differentiation across the sample, which we relate to the evolutionary 
state of the cloud. We relate our observations to the cloud physical, kinematical and evolutionary 
properties, and also compare them to the predictions of simple chemical models. The implications of 
this work for understanding the origin of the clumpy structures and chemical differentiation observed in 
dense clouds are discussed. 



1 Introduction 

Low-mass stars form in dense cloud cores. Understanding how stars form thus requires detailed knowledge 
of the physical and chemical evolution of molecular clouds. The classic picture, in which star formation is 
magnetically-dominated and mediated by ambipolar diffusion [e.g. 1 , posits long molecular cloud lifetimes 
(~ 10 7 years). However, there is an increasing realisation that dissipation of turbulence, driven externally 
at large scales, plays a vital role [e.g. H]. In this case, star formation is dynamic and rapid, and molecular 
cloud lifetimes are short (~ 10 5 — 10 6 years). Dense cloud chemistry will follow different evolutionary paths 
depending upon the lifetime of molecular clouds, and hence on the time-scale of low-mass star formation 
[e.g. 3 . One may also expect that the spatial distribution of molecules should also reflect the underlying 
cloud dynamics. Observations of dense interstellar clouds in a variety of molecular tracers show that they 
contain a distribution of dense gaseous structures [I] . These structures exhibit a size distribution ranging 
from that of dense cores [~ 0.05 — 0.1 pc, 5 , to clumps [~ 0.01 — 0.05 pc, O [7] down to so-called small 
clumpinos [~ 0.005 pc, H]. The spatial distributions of molecules in cold, apparently quiescent, clouds do 
appear to show chemical abundance gradients, i.e. chemical differentiation. Previously, the best studied 
sources have been the clouds TMC-1 and L134N [21 [TU1 HI] • These clouds show striking differences in the 
location of the emission peaks in several molecules. In TMC-1, for example, emission from the carbon 
chain molecules (e.g. the cyanopolyynes) is observed to be anticorrelated with that of ammonia and 
N2H+ and several deuterated species [e.g. [T2] [T31 [H]. These clouds also show distinct emission peaks 
from other specific molecules, such as methanol and SO. High resolution observations show that chemical 
differentiation is present at all spatial scales [15j [16]. Although, in dense cores, depletion of CO, CS 
and other heavy molecules on to dust grains is undoubtedly important |17l 118] , the reason for the spatial 
differentiation seen amongst other molecules is unclear, as is its prevalence in other molecular clouds. 

It is possible that there is a link between the observed chemical differentiation and the clumpy struc- 
ture, and hence to the physical origin of clumps through the cloud turbulence. There are currently two 
possibilities [see 19, 20 . First, numerical simulations show that dissipation of externally driven turbu- 
lence can produce structures, similar to the clumps and cores that are observed, through the process of 
'turbulent fragmentation' [2J. Second, winds and outflows have a significant effect on the structure and 



Table 1: Core Properties 
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evolution of the surrounding protostellar envelope, containing energies sufficient to physically disrupt the 
envelope material and to alter the chemistry. These outflows begin very early in the star formation pro- 
cess, and cloud turbulence could be excited by them [2T] . In this case, the cloud structure is regulated by 
the action of stars [e.g. [22]. The presence of parsec-scale flows [e.g. [25] and wind blown bubbles [21] are 
clear observational evidence that protostars affect their environment in this way, and that the influence of 
existing star formation can have a very long reach. In this picture, the turbulence should be a function of 
the local star-formation activity, however it is unclear if enough kinetic energy can be injected in all clouds 
to sustain the observed level of turbulence. As the nature of the turbulence should affect the formation 
and evolution of clumpy structure, we might therefore expect that there should be chemical differences 
between starless cores and those containing protostars at the Class and/or Class I stages [25], with the 
former being the product solely of turbulent fragmentation. Models have previously been developed of 
the spatial and temporal chemical evolution in clumpy molecular clouds [e.g. [IT] [25] [27] [28, 29, 30] I3T] . 
but with limited success. 

To quantify the chemical effects that embedded protostars may have on cloud structure and chemistry, 
we have previously mapped the Barnard 5 cloud (B5) in many molecular lines. B5 is an active region of 
star formation in Perseus and contains molecular outflows, wind-blown bubbles, and regions where the 
outflows appear to be interacting with dense clumps [52] [33] [33] . Most previous observations have focused 
on the region containing the Class I source, B5 IRS1, and the molecular distributions over the rest of the 
cloud were previously poorly known. These observations revealed a high degree of spatial differentiation, 
with similarities to the molecular distributions in TMC-1 [35], but also to some other starless cores in 
Taurus [36]. Comparison with dynamical-chemical models of B5 [26], based on the Norman & Silk [22] 
picture, appear to indicate that molecular ices are continually being formed and destroyed [7J 136] . but 
that processes not included in these models must be responsible. However, B5 is only one such source; 
to constrain any chemical models multi-molecule maps of more sources are needed, for example, to see 
if similar molecular differentiation occurs in dense clouds with varying degrees of star formation activity 
and protostellar evolution. We have therefore increased our sample to observe dense core chemistry in 
different natal environments: starless cores, and those containing young stellar objects (YSOs) at the 
Class and Class I phases. The sources were selected from the sample mapped by Caselli et al. [37] 
in N2H" 1 " at 0.063 km s _1 , 54" resolution; lower resolution ammonia maps also exist [5J. Caselli et al. 
[57] found that the position of the N2H+ peak emission tended to be close to the position of peak NH3 
emission but generally not coincident with that of the star. From our B5 observations, we found that 
several molecules are strongly anticorrelated with N2H 4 " and NH3 emission at small scales, as well as 
with each other (e.g. methanol and the carbon chains). Hence, the N 2 H + maps of Caselli et al. [37] can 
be useful to predict where other molecules could be expected to peak in a cloud, assuming the chemical 
differentiation has the same underlying cause. 

The six cores we observed have a range of star-forming activity (see Table [l}. TMC-1C is a starless 
core, with a region of diffuse sub-mm emission to the south, which extends in a ridge beyond our map to 
the north [38] . This cloud has previously been mapped in carbon chain molecules |39] and in methanol 
and HC 13 0+ [5]. L1389 is associated with an IRAS source, IRAS 04005+5647, although it has been 
classified as pre-stellar, or an extremely young protostellar source 40] . L1262 contains both a stellar 
core, IRAS 23238+7401, which drives a protostellar outflow, and a starless core [41 . L1251E contains 
a YSO, IRAS 22385+7457, which drives an outflow, as well as two infrared sources likely to be T Tauri 
stars 42] . There are also several infrared sources in the immediate vicinity of our map. L1512 is a starless 
core, containing a north-south ridge of diffuse sub-mm emission |41j . Per7 contains a sub-mm peak in 
addition to the Class I source, IRAS 03295+3050 [23]. 
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Table 2: Molecular Line Properties 
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We chose to observe four molecular species (see Table 2). The isotope C 18 serves to trace much 
of the low density (~ 10 3 cm -3 ) material in the cloud. Methanol is formed solely on grains in dense 
clouds [e.g. 33]. Once desorbed into the gas, methanol molecules only have a finite lifetime until they are 
destroyed. Hence, the CH3OH distribution provides an excellent observational measure of the degree of 
gas-grain cycling in clumpy clouds. Cores showing weak methanol emission will tend to be older since 
any CH 3 OH liberated previously will gradually re-accrete back onto the dust. As we discovered a new 
'cyanopolyyne-peak' position in B5, similar to that in TMC-1 |35] , we also mapped the cyanoacetylene 
emission in these clouds. The formation and destruction of cyanopolyynes can provide information on 
those molecules returned to the gas, and also define chemical time-scales in dynamically-evolving clouds 
[29] . Markwick et al. [29] concluded that such chemically interesting positions should appear to be offset 
from the position of any YSO present; as well as from the peak ammonia emission. The detection of these 
organic molecule 'hotspots' is thus only possible by making large-scale HC3N maps, to include positions 
distinct from the nominal (0,0) positions of previous searches [cf. [57] ■ The HC3N transition chosen 
traces gas at densities above ~ 10 5 cm -3 . Finally, to explain the large number of S-bearing molecules 
detected in dense clouds, chemical models infer that elemental sulphur is only modestly depleted in dense 
clouds, relative to more refractory elements [45] . Interstellar sulphur chemistry involves primarily neutral 
reactions, and so it defines yet another, albeit slow, chemical time-scale. Shock chemistry or grain-surface 
catalysis have been proposed to be important for S chemistry [30 , 46, ; the latter mechanism is supported 
by the detection of OCS in interstellar ices [UJ. We therefore mapped the Jk — ^2 — 2i line of SO. 

The plan of the paper is as follows. In Section 2 we describe the observations and how the data 
analysis was performed. In Section 3, we detail the chemical morphology of each source and note the 
general trends we see in chemical differentiation. Simple chemical models of the chemical evolution in 
clumpy media are presented and compared to the observational data in Section 4. Conclusions from this 
study are given in Section 5. 

2 Observations, Data Reduction and Analysis 

Single dish maps of HC3N, SO, CH3OH and C ls O transitions at 3mm were made using the Onsala 20m 
telescope, Sweden, in April 2005. With the exception of CH3OH, the spectra were observed in frequency- 
switched mode, using 12.5 kHz channels to provide a velocity resolution of ~0.04 km s -1 . CH3OH was 
observed with the same resolution, but in position-switched mode due to the number of transitions in 
the band. The 20m beam size at these frequencies is approximately 40", and the main beam efficiency, 
r\ m b, is 0.46. Table [2] lists molecular line details for the transitions that we observed. 

The cores observed were taken from Caselli et al. [37] , and include cores at a range of evolutionary 
stages, from prestellar to Class I (Table [TJ. The areas mapped cover the N2H + emission regions within 
each core [37] , ranging in size from 3-20 sq. arcmin. The typical RMS noise level in the maps is 0.1-0.3 K. 
Contour plots and channel maps are shown in Figs |3.1| |3.1[ 

Data reduction was carried out using the radio data reduction packages XS (bespoke software for the 
Onsala 20m telescope), CLASS and SPECX. Single Gaussian fits provide a good fit to the line profiles. 
For CH3OH, the transitions detected in each spectra were fit with dependent Gaussians. The velocity 
and linewidth of each transition was fixed with respect to the 2o ^lo(^4 + ) transition at 96.741 GHz, since 
the separation of the K transitions is fixed by the structure of the molecule. The results of the fitting 
procedures are given in Table § and for CH 3 OH, values are given for the 2 — >lo(A + ) transition. For 
lines that were not detected, the 3cr error estimate is given as the upper limit to the integrated intensity. 
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Since we have observed single transitions in C 18 0, SO and HC3N, and detect only one or two CH 3 OH 
transitions in a large fraction of the spectra, we determine a lower limit to the column density |48j . In 
this procedure, the rotation temperature can be approximated as T rot = for non-linear molecules 
(CH3OH), and as T rot = ^ for linear molecules (C 18 0, SO and HC3N). A lower limit to the column 
density can then be calculated as: 

^ /T^dvQ^e 3 / 2 non- linear (1) 
^ /T mb dvQ(^)c linear (2) 

where T m b has been calculated from using r\ m h-, partition functions (Q[T]) and statistical weights (g u ) 
have been taken from the JPL molecular line database [39], energy levels and Einstein A values (A„;) 
have been taken from the Leiden Atomic and Molecular Database j5D] and Cragg et al. [ST], listed in 
Table[| For CH 3 OH, the partition function Q(T rot )=1.28T^ was adopted from Menten et al. [52] , and 
takes into account equal populations of A and E species |53j . 

For sources where we have detected 3 or 4 CH 3 OH transitions, we have also carried out a rotation 
diagram analysis. The resulting column densities generally agree to within a factor of a few. However, 
as we have detected only one strong transitions in each A and E species towards most positions, we have 
based our CH3OH analysis in this paper upon the 2q — > lo(A + ) transition, in the same manner as for the 
other molecules. This provides a coherent analysis of this data that we can compare with simple chemical 
models. 

From our C 18 observations, we can then derive abundances relative to H 2 following Frerking et al. 
[M] for N(C 18 0) > 3 xlO 14 cm" 2 : 

x 10 21 cm" 2 (3) 

We derive H2 column densities ~ 2 x 10 22 cm~ 2 within the clumpy regions of our data, leading to 
abundances on the order of x(CH 3 OH) - 2 x 10" 8 , x(HC 3 N) - 8 x 1CT 11 and sc(SO) - 4 x 1(T 10 in 
the molecular clumps, or 'hotspots'. 
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3 Molecular Morphology 

We now describe the observed chemical morphology. We first discuss the trends between the four 
molecules in each individual source, focusing on the emission present in the position- velocity channel 



maps (Fig. 3.1 -3.1 1. We then summarise the apparent general trends. 



3.1 Individual Sources 
L1512 

L1512 shows several extended clumps and ridges of emission in C 18 0. The southern part of the map 
shows two clumps ((-60",-30") and (0",-30")) which surround the compact HC 3 N clump (-30",-30"). 
Emission from SO is also seen in extended ridges and clumps. The strongest clump is associated with 
the single HC3N clump (-30", -30"), while the second strong clump is associated with the CH3OH peak 
to the north (-60", 120"). 

We see SO emission at 6.62 km s _1 with none apparent in either CH3OH or HC3N. The strongest 
SO emission clump appear in gas in which CH 3 OH is present, at 6.88 km s _1 , the SO and CH 3 OH 
emission peaks are offset by approximately 30". The methanol peak is at 7.12 km s _1 , where most of 
the weaker clumpy emission is found; at this velocity, there is no corresponding spatial SO emission in 
this region. There are a couple of weak CH3OH emission clumps present to the south, and these are 
embedded in the same gas as two SO clumps. An HC3N clump shows up at -6.88 km s _1 that is spatially 
distinct from the SO and CH 3 OH clumps but partially overlaps the weaker extended SO emission. The 
main HC3N clump is at 7.12 km s _1 , overlaps the extended SO and CH3OH there, but is distinct from 
the weak SO and CH3OH clumps, and is located about 100" from the major SO-CH3OH clumps. In 
fact, we find that both these HC3N clumps coincide with minima in the C 18 distributions, as well as 
with the N 2 H + peak (Vlsr = 5.27 km s" 1 ), perhaps indicating that molecular depletion is important 
at these positions. However, the fact that a methanol clump is present at 7.12 km s _1 in a region with 
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Figure 1: L1512. Top: Contour maps of integrated intensity, from left: C 18 0, SO, HC 3 N and CH 3 OH. 
Contours are 50%-100% in 10% intervals of the peak integrated intensity in K km s _1 , and la rms is 
~ 0.04 K km s™ 1 . A cross marks the N2H" 1 " peak [37]. Bottom: Channel maps of integrated intensity. 
Contours levels are marked on each plot. Clockwise from top left: C 18 0, SO, CH 3 OH and HC3N. 

little/no C 18 emission somewhat confuses this interpretation. C 18 emission peaks are anti-correlated 
with emission peaks in all other molecules. 

In summary, we see gas with only SO, gas with strong SO and CH3OH emission intermingled but with 
distinct emission peaks (i.e. clumps), as well as gas with strong CH3OH emission and no/weak SO. The 
HC3N clumps are markedly anticorrelated with SO and CH3OH clumps in position and are correlated 
with C 18 troughs. These could be cyanopolyyne late-time depletion peaks [55]. 

TMC-1C 

TMC-1C shows several clumps of emission in C 18 extending in a north-south ridge. The peak of emission 
is near the centre of the diffuse sub- mm emission (-30", -60"), with several smaller clumps leading north 
to a more extended emission region at the edge of our map (-90", 120"). HC3N emission is confined to 
the southern half of our map, where it appears in a ring, part of which overlaps the C 18 clump. SO 
emission appears in a relatively narrow ridge of emission extending to the north-west. Emission from 
CH3OH has a similar morphology to emission from SO, but is more extended, and peaks further to the 
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Figure 2: TMC-1C. As Fig. |3.1| A triangle marks the approximate centre of the diffuse sub-mm emission 
to the south, which extends slightly to the east, and to the north beyond our map [38]. A cross marks 
the N 2 H+ peak [37]. 



north. 

In this source SO displays very clumpy emission at 5.f2 and 5.38 km s _1 ; there is only weak/no 
CH 3 OH emission at 5.f 2 km s _1 , where most of the SO molecular clumps are. At 5.38 km s _1 , CH 3 OH 
has 3 clumps and there is evidence for coincident weak SO emission. There is an SO clump coincident 
with a minimum in C 18 emission; this is approximately the location of N2H + emission maximum |37j . 
and suggests a depletion zone here. This is curious: in L1512 we have HC3N correlated with a possible 
depletion zone, whereas here it is SO and there is only weak/no HC3N present - or at least tending to 
avoid it. The N 2 H + peak is at Vlsr °f 5-27 km s _1 , at a similar velocity to emission from C 18 0, CH 3 OH 
and HC3N. The velocity of SO emission at this position is blue-shifted by 0.1 km s _1 , which presumably 
means that the depleted material does not contain the SO peak. The structure of the C ls O, HC3N and 
SO emission could be interpreted as that of a clumpy, fragmented torus in 3D. In this scenario, the SO 
emission is actually coming from a chemically differentiated fragment of the torus/ring. The HC3N peaks 
are spatially distinct from the CH 3 OH clumps at 5.38 km s _1 but there is weak SO emission in them; an 
exception is the SO clump near (0,+60") which appears to have no HC3N. Again, C 18 emission peaks 
are anti-corrclatcd with emission peaks in all other molecules. 

Thus, in the SO clumps, methanol is cither deficient or absent. To some degree, SO and HC3N 
appear to be intermingled, although the clumps/peaks are physically distinct, and one SO clump has 
no HC3N. The CH3OH clumps probably have some SO in them but are spatially distinct from the 
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L1262 

Emission from C 18 towards L1262 is seen in two strong clumps, one just to the east of of the IRAS 
source (0",0"), and one just to the east of the sub-mm source (60",-30"). Emission from HC3N is seen in 
a clumped ridge, with one clump, like that of C 18 0, to the east of the sub-mm peak, and one associated 
with the IRAS source, west of the C 18 clump (30",-30"). SO emission peaks in a clump to the north 
of the sub-mm peak (0",30"), and extends in a clumpy ridge to the north and south. A further clump is 
seen to the west, at the edge of our map, where the C 18 emission also extends. Emission from CH3OH 
is seen in an extended clump associated with the sub-mm source (0",30"), and also in a compact clump 
to the east of the IRAS source (60", -30"). The main SO and CH3OH clumps are co-incident with the 
HC3N clump that is near the sub-mm peak. 

There is an SO clump at 3.62 km s _1 that appears to grow and merge into the major SO clump at 
3.88 km s _1 , where CH3OH emission becomes more evident. At 4.12 km s _1 , the SO emission clearly 
shows up as ~ 4 clumps. There are 2 CH 3 OH clumps and these are evident at 4.12 and 4.38 km s _1 . In 
both cases, the CH3OH appears to be mixed with the enhanced SO emission, although at 4.38 km s _1 
the SO emission is absent. The strongest SO and CH3OH emission are to the north - the clumps/peaks 
are spatially distinct and should have different SO/CH3OH ratios. The methanol clump at 4.38 km s _1 
is almost 'pure' CH3OH - emission from other molecules is weak/absent. There are 2 HC3N clumps with 
the strongest to the north at 3.88 km s _1 . This clump is associated with the strongest SO peak. Only 
the weaker HC3N clump to the south is also evident at 4.12 km s _1 which is located across a region of 
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weak/absent SO and CH3OH emission. There are two C 18 clumps, which are distinct both spatially 
and kinematically. Through both of these, we see a kinematic sequence of C 18 and SO, then C 18 and 
SO plus CH 3 OH, then just CH 3 OH as we move from the blue to the red in velocity channels, although 
the spatial peak in emission from these molecules are offset from each other. 

In LI 262, we detected 4 SO clumps. The strongest one also contains quite strong HC3N and weak 
CH3OH emission at 3.88 km s _1 ; however, at 4.12 km s _1 the CH3OH clump grows stronger and more 
distinct while the HC3N emission disappears. This is similar to that seen previously, where the strongest 
SO and CH 3 OH clumps tend to be in gas where the molecules coexist. We note that the methanol 
clumps are spatially anticorrelated with the positions of the HC3N ones, and that the weaker CH3OH 
and HC3N clumps have little SO. However, the positional offset of the emission peaks of each molecule 
is only marginal (~0.5 beamwidth). 

The SO, CH3OH and HC3N peaks are all anticorrelated with C 18 peaks. These molecules also tend 
to avoid the region around the N 2 H + emission peak (at Vlsr = 4.11 km s _1 ), suggesting that molecular 
depletion on grains could be influencing the observed morphology. There may be a kinematical sequence 
through the clumpy structures around (-60, 0). 

Per 7 

Per7 shows several regions of clumpy, asymmetric emission from all four molecules, most of it contained 
in ridge extending to the north-east and south-west from the sub-mm peak. Emission from C 18 is 
relatively bright towards this source, while emission from HC3N, CH3OH and SO is fairly weak. C 18 
peaks in an extended region to the NE of the sub-mm peak (around (30", 30")). HC3N peaks in the 
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same place, but emission is much more compact, with a second emission peak to the SW (-90", -60"). SO 
appears in a series of clumps along the NE-SW ridge, with the brighter clumps those to the SW (-90",- 
60"). CH 3 OH also appears in a series of clumps along the NE-SW ridge, but the brightest clump is to 
the far NE (90", 90"), where no SO emission is seen. Another peak occurs in between the two strongest 
SO emission clumps (-60", 0"). 

Very clumpy structure is evident in emission from SO. At 6.38 km s _1 there are 2 SO clumps to the 
north with no/little CH3OH or HC3N emission. However, the peak emission clumps of SO and CH3OH 
are in the same gas, but physically distinct (there is also a possible/weak HC3N clump at this velocity). 
At 6.62 km s _1 , SO and CH 3 OH show several clumps in gas where they are cominglcd but with the peaks 
positionally offset from each other. There are 3 major CH3OH clumps and it seems that the northern 
one contains little or no SO and definitely no HC3N. At 6.62 km s _1 , there is one clear HC3N clump that 
actually sits in the CH 3 OH minimum. In fact, once again, the HC 3 N and CH 3 OH emission is generally 
anticorrelated. The C 18 does not show compact peaks in the contour map, but is also anticorrelated 
with HC3N and SO. All the molecules avoid the associated N 2 H + peak. 

In this source we again see clumps with essentially only SO present, as well as that the peak emis- 
sion from SO and CH3OH emanates from gas in which they are comingled, but where their emission 
peaks/clumps are spatially distinct. There appear to be methanol clumps with no HC3N and little or no 
SO. Again, there appears to be a chemical gradient through the clumpy structures at around (+100, 0): 
6.38 km s" 1 : SO only - 6.62 km s" 1 : strong CH 3 OH and unchanged SO : 5.88 km s" 1 strong CH 3 OH 
no SO. 

L1389 

L1389 shows relatively compact and weak emission from all four molecules. C 18 peaks in a clump to 
the east of the IRAS source (-30", 0"). HC3N peaks in a clump south of the C 18 emission. Emission 
from SO also peaks to the east of the IRAS source (-30", -30"), but is more extended. CH 3 OH emission 
appears in a strong clump to the south of the IRAS source, extending to the west (-60", -30"). 

Chemically, this is a morphologically simple cloud. The SO shows a single clump at -4.88 km s _1 
with weak, diffuse emission at -4.62 km s _1 ; this clump appears to have very weak HC 3 N emission but 
no CH3OH. The CH3OH clump is at -4.62 km s _1 and this corresponds to a region of relatively weaker 
SO emission. The HC3N clump appears at -4.62 km s _1 and and is spatially distinct from the SO peak 
in velocity and the CH3OH peak in position. 

So, there is a clump in which SO is prevalent, although it does contain some HC3N. A distinct 
HC3N-CH3OH anticorrelation is also evident. All the molecular clumps (positions of peak emission) are 
located away from N 2 H+ peak, supporting the idea that molecular depletion is important here. 

L1251E 

Towards L1251E, emission from C 18 is seen in several small clumps, the strongest of which is associated 
with the IRAS source (-210", 0"). Two clumps are also seen to the north and south of an infrared source, 
either side of a clump of HC3N emission (-150", 0"). HC3N emission is seen in a second clump associated 
with the IRAS source (-210", 0"). SO emission is seen in a ridge of emission extending to the north 
(-120", 120"), not clearly associated with emission from any of the other molecules, nor any of the infra- 
rcd/IRAS sources. The clumps near the second infra-red source are associated with the weaker CH3OH 
clumps. CH3OH emission is seen strongly peaked in a clump associated with both the IRAS source and 
the first infrared source (-180", 0"), with a second peak to the north, just south of the main SO clump 
(-90",60"). 

The strong SO emission clump to the north, evident from -4.88 to 4.12 km s _1 , contains little or no 
CH3OH, no HC3N, and is unconnected to the level of C 18 emission present. The CH 3 OH emission is 
extended and fragmented and three main clumps are evident. There are 2 HC3N clumps present. The 
HC3N peak/clump at (-120, 0) appears quite extended, shows weak CH3OH emission but none from SO. 

In summary, there is an SO clump that contains a small amount of HC 3 N emission and only a small 
amount of CH3OH emission. The HC3N - CH3OH peaks are anti-correlated to the east, where the 
HC3N peak appears in a region with no CH 3 OH emission. To the west, HC3N and CH 3 OH emission are 
comingled, but the peaks are offset. 
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Figure 5: L1389. As Fig. 3.1 A star marks the IRAS source, a triangle marks the sub-mm peak [40 
and a cross marks the N 2 H + peak [57] , 



3.2 General Trends - Summary 

Several trends in the chemical differentiation are evident in our maps. The morphology throughout 
our sample ranges from one or two isolated clumps detected in each molecules, to extended emission 
containing a strongly-emitting clump and a several weaker ones. 

SO and CH 3 OH Clumps 

We find clumps of gas that contain significant emission only from SO (e.g. L1389, L1251E, PER 7, 
TMC-1C), with only one of these showing significant emission from HC3N (L1262). There is evidence 
that some SO emission maxima could correspond to depleted regions (TMC-1C). Strong extended SO 
emission is also evident where CH3OH emission becomes stronger; in this case we find that the dominant 
CH3OH and SO clumps can exist in the same gas but that they are always offset from each other. On 
the other hand, there are some methanol clumps from which SO emission appears to be absent. Based 
on comparing maps in velocity space, there is some kinematical evidence (e.g. in L1262, PER 7, L1512) 
for well-defined emission sequence in which there are : regions with only SO, then a region with SO 
and CH3OH intermingled - then a region with only CH3OH, suggesting we could be looking down a 
cylinder/ filament . 
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Figure 6: L1251E. As Fig. 3.1 The IRAS source is marked with a filled star, and unidentified IR sources 
with empty stars [12J. A cross marks the N2H+ peak [3"T] . 



The HC 3 N Clumps 

We find that the HC3N clumps are always spatially anticorrelated with the CH 3 OH clumps. Some of 
these clumps can contain some very weak CH3OH emission, and moderately strong SO emission (L1262), 
but there are others where the HC3N is the sole molecule present (L1251E). 



Depletion: C 18 O Distribution & A^-ff -1 " Cores 



C O emission is tracing the less dense gas towards these sources, since the critical density of the transition 
we have observed is 1 or 2 orders of magnitude less than the transitions we observed in SO, HC3N and 
CH3OH. The emission generally has several clumps embedded within an extended, diffuse region, and 
emission from spatially distinct regions is generally seen to be kinematically distinct as well. For example, 
towards L1512, emission to the north-west is seen at 6.88 km s _1 , while emission to the south-east is seen 
at 7.12 km s _1 . There is a clear anti-correlation between emission from C 18 and all other molecules 
towards the youngest sources, those classified as prestellar cores. Towards the older sources, where 
protostars have formed, emission becomes more intermingled, and C 18 emission can be seen associated 
with SO and CH3OH peaks, although the spatial positions are offset (e.g. L1262). 

Another population of objects implicitly present in our maps are the N2H + cores mapped by Caselli 
et al. [37 . These are probably the most stable, longest lived cores in the clouds, apart from the ones 
that apparently harbour protostars, and are almost certainly sites of extensive molecular depletion onto 
grains. All molecules in our sample tend to avoid these positions. An exception is in L1512 where the 
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HC3N clump is coincident with the N 2 H + peak in both position and velocity. 



4 Chemical Differentiation in a Clumpy Medium: Theory 

We conclude that we are observing chemical morphologies which are produced by a time-dependent 
gas-grain chemistry which is evolving in a turbulent dense medium. Clumps form and dissipate, perhaps 
coalescing, and long-lived stable structures (i.e. dense cores) are produced in which star formation occurs. 
In this section we show how simple chemical models can help explain these observations. 

Clumps are formed out of the intcrclump medium (nn ~ 10 3 cm -3 ), perhaps by compression in 
shock waves or in MHD waves [SI]. Our observations indicate that molecular desorption/deplction 
from/onto dust grains plays an important role in determining the chemical differentiation. For a molecule, 
X say, of molecular weight fix, the depletion (accretion) time-scale at 10K is 

1/2 

t x = 2.2 x 10 9 -^- yr (4) 
n H 

where a total grain surface area per H nucleon of 3 x 10~ 22 cm 2 and unit sticking efficiency have been 
assumed [cf. [57]. In these dense regions (nn ~ 10 4 — 10 5 cm~ 3 ) accretion and chemical reactions are 
driven more rapidly and dynamical fate of the clump becomes an issue. It is therefore useful to know the 
relevant physical time-scales that arc important for clump evolution. 



4.1 Clump Physics 

The fate of clumps can be related simply to parameters that can be derived from observations. From our 
observational data we can obtain estimates for the clump mass, M, its virial mass, M v ; r , and radius R. 
With R known, the virial mass for a sphere of constant density can be estimated from [SS] 

M vir ~ 210(Az/) 2 v i? M G (5) 

where (Aw) av is the average line-width of the clump gas (H2 and He) in km s _1 and R is in parsec 
For clumps in molecular clouds, there are three possible evolutionary tracks [cf. [6]. A clump will be 
stable against gravitational collapse if (0.5M vir < M < M vir ), whereas if it is gravitationally unstable 
(M > M v j r ) it will collapse on the free-fall time-scale 

4.35 x 10 7 

ts = — 172 — y f ( 6 ) 

n H 

If (M < 0.5-Mvir) then the clump is gravitationally unbound and so it will dissipate on a time-scale 

2R 

£diss ~ 7^— (7) 

where C e g is the effective sound speed in the gas. The fate of some clumps can also depend upon the 
local environment. A cloud volume containing a number of clumps of filling factor 7, and line-of-sight 
relative velocity v rc \ can collide and coalesce (i.e. merge) with a time-scale 

2R f 1 \ 

Emerge ~ I ] 75 — 1 I (°) 

Clumps with any dynamical fate can grow when the merger time-scale is sufficiently short; unbound 
clumps can do so when this is comparable or less than the dissipation time-scale, i.e. when t merge < idiss 
[3 US]- If we take the derived clump properties of the Taurus clouds (TMC-1C and L1512) as an 
example, we then, following Peng et al. [6], we estimate that typically R ~ 0.02 — 0.03 pc and the clump 
masses are in the range M ~ 0.3 — 1 M Q . A detailed analysis of clump properties is beyond the scope of 
this paper, however, we note that throughout our sample we can identify clumps that are unbound, as 
well as some that are gravitationally stable and others that are unstable. Some clumps, such as CH3OH 
and SO clumps in TMC-1C and L1512, may actually show substructure at higher spatial resolution. In 
this case, such smaller clumps would be more likely to be unbound and, given their proximity, perhaps 
more prone to interaction. 
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4.2 Chemical Model 



We now show how a simple static model of clump chemical evolution, with accretion on to dust and 
selective injection of simple mantle molecules, can explain the observations. 

As discussed in section 1, methanol cannot be produced efficiently in the gas phase, and is thought to 
be formed on the surfaces of dust grains. Hence, the CH 3 OH-rich clumps in our sources must necessarily 
trace clumps in which grain ice mantles have recently been liberated into the gas phase. The exact 
mechanism that causes the desorption of the ices is unknown but we expect that it will be connected to 
clump evolution. Markwick et al. |29j proposed that the desorption is driven by grain-grain collisions, 
induced by MHD waves generated by clump motions. Alternatively, Dickens et al. [65] speculated that 
grain-grain streaming resulting from clump collisions may be responsible. In this paper, we make no 
assumptions about the dynamical mechanism, other than that it is related to clump formation, and so 
ice mantle liberation, if it occurs, happens at the same time that the clump is formed. This is obviously 
a simplification since there must be a period of CO accretion to provide the precursor ices for CH3OH 
production. We assume that a clump forms instantaneously from the interclump medium, and simply 
follow the chemistry until all the molecules arc condensed on grains. In reality, the chemical evolution 
will be sensitive to the dynamical evolution of the clump (e.g. whether it is collapsing, coalescing or 
dissipating). We do not consider the clump physical evolution explicitly in this static model; this is 
discussed in detail elsewhere [55] , 

Based on the critical densities of the observed molecular tracers, we adopt a value of jih = 2 x 10 5 cm~ 3 
as representative of the clump densities, which we consider to be approximately spatially constant. This 
choice yields tg ~ 1 x 10 5 years, and tdiss ~ 2 x 10 5 years, and clump masses similar to those derived 
from the observations, viz 

M = O.35M ( — ( - — \ (9) 

vo.02 pcy v 2 x 10 cm / 

We assume a cosmic ray ionization rate of 3 x 1CP 17 s _1 , a temperature of 10 K, and a visual extinction 
of 10 magnitudes, which is sufficiently large to render photo-processes effectively irrelevant. Freeze-out 
of molecules onto grains occurs at a rate given by eqn. Q, except for N2 which is assumed to remain in 
the gas phase [e.g. [T7]. The chemical model is based on the model of Rodgers & Charnley [50], and has 
been updated to incorporate recent laboratory measurements of the dissociative recombination channels 
for several ions, including CH 3 OH+ (44j EU E21 E31 HE] • 

The initial molecular abundances in the newly-formed clump are determined by those in the interclump 
medium, which are themselves controlled by both (i) the chemistry in these regions, and (ii) whether 
the material previously passed through a dense phase. This latter effect may be important, since if 
most clumps eventually dissipate, recycling dense cloud material back into the surrounding gas, the 
abundances in the interclump medium retain a memory of the higher density phase [cf. [26] . We assume 
that the material in the interclump medium is predominantly molecular, and that « 99% of the carbon 
and nitrogen incorporated into a clump is initially present as CO and N 2 , with w 1% in neutral atomic 
form. The oxygen not bound up in CO is assumed to be atomic, as is all the initial sulphur. To account 
for depletion, the gas-phase sulphur abundance is assumed to be 10 -8 |45j . and we assume complete 
depletion of all metals. The initial fractional ionization is 1.5 x 10~ 7 . 

4.3 Results 

Based on the relative binding energies of species thought to be present on interstellar grains [67], we 
consider three types of clumps: 

• Type I. No sublimation of any ice species. The initial gas-phase composition is equal to that of the 
interclump medium. 

• Type II. Sublimation of the most volatile ice species. We assume injection of CO2 and H2S with 
abundances of 10~ 5 and 10 -7 respectively. 

• Type III. Sublimation of tightly-bound species. We inject CO2, CH3OH, and H2O with abundances 
of 10~ 5 , 2 x 10~ 8 , and 3 x 10~ 5 respectively. 
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Figure 7: Abundances versus time for selected species. The vertical dotted lines mark the time-scales for 
gravitational collapse (tg) and clump dissipation (tdiss)- Top: type I cores (no ice sublimation). Middle: 
type II cores (CO2 and H 2 S sublimation). Bottom: type III cores (CO2, H 2 0, and CH 3 OH sublimation). 



The abundances of sublimated CO2 and H 2 are based on observations of interstellar ices in dark clouds 
[55] . and the CH3OH abundance is chosen to match our derived column densities (see section 2). H 2 S has 
not been detected in interstellar ices, but based on the fact that many of the most abundant ice species 
are simple hydrides, it is likely to be present; our assumed abundance is consistent with the observed 
upper limits [68] . Hydrogen sulphide has a low sublimation enthalpy [67] , and so it will be released along 
with C0 2 in type II cores. We neglect H 2 S injection in type III cores in order to distinguish between the 
amount of SO that can be formed from H2S, and the amount that can be formed by simply boosting the 
OH abundance (as will occur when water is injected). 

Figure |4~3] shows the time evolution of each of the four observed species - CO, SO, HC3N, and CH3OH 
- in each clump. Also shown are N 2 H + and S02- The free-fall time, tg, and the clump dissipation time, 
fdiss are marked on each panel. In all cases, molecules are depleted from the gas on a freeze-out timescale 
of a few xlO 5 yr. For type I clumps, SO and HC3N are produced by chemical reactions in the dense 
gas, peaking between 5 x 10 4 and 10 5 yr. SO is formed via the reaction of atomic S with OH, and 
HC3N from CN + C2H2. Type I clumps younger than ~ 10 5 yr can account for regions with both SO 
and HC3N emission but no methanol. At longer times, SO is destroyed faster than HC3N, so long-lived 
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clumps of this type are able to explain those regions which are HCsN-rich but SO-poor. In type II 
clumps, the injection of H^S significantly enhances the amount of SO produced in the gas. Compared 
to type I clumps, the HC3N peak is slightly reduced, and shifted to later times. This is due to the 
reduced hydrocarbon abundances which result from the fact that the injected CO2 drives down the C + 
abundance. Type II cores can explain regions with large SO abundances, together with small HC3N 
abundances and no CH 3 OH. Type III clumps are obviously the only ones which can explain CH 3 OH-rich 
clumps, but we find that the methanol only survives for ~ 10 4 yr before it is destroyed by protonation 
followed by dissociative recombination. The HC3N in these clumps only peaks after all the CH3OH has 
been lost, which naturally explains the anti-correlation between methanol and cyanoacetylene. The SO 
abundances in type III clumps are intermediate between those of types I and II, and peaks occur at 
around 2 x 10 4 yr. At this time, much of the initial CH3OH has been destroyed, and no significant 
amounts of HC3N have been formed. Therefore, type III cores younger than ~ 10 4 yr will contain both 
CH3OH and SO, whereas slightly older cores will only contain SO. Because the injected water leads to 
an increased OH abundance in the gas, type III cores eventually evolve to a state with SO2/SO > 1, so 
the SO/SO2 ratio may potentially be used to place more firm constraints on the ages of these clumps 

[&f. ESI ESI. 

4.4 Application to Observations 

We can compare our model results with the observational dataset, to try to place constraints on the 
ages and the dynamical states of specific objects. Considering the general trends described in section 
3.2, we find that our model naturally accounts for the CH3OH-HC3N anti-correlation. Clumps with 
CH3OH emission must have ages of t< 10 yr, and those with no SO emission must be even younger 
(t< 5 x 10 3 yr). Because these ages are much less than those associated with the dynamical evolution 
of the clumps, the dynamical state of methanol dumps cannot be constrained from the observations. As 
these type III clumps evolve, they eventually lose all the initial methanol; for 3 x 10 4 <, t < 10 5 yr, they will 
appear only in SO, after which time the HC3N abundance rises. The evolution of type III clumps provides 
a simple explanation for the observed transition from pure-CHsOH to mixed-CH30H-SO to pure-SO gas 
in different velocity channels. Clumps with both SO and HC3N emission could be old type II or III 
(t > 10 yr), or possibly younger type I clumps. Regions containing only HC3N emission are likely to be 
old type I clumps, with ages f~2x 10 5 yr ~ idiss- As these latter clumps have ages comparable to tdiss; 
it is probable that pure HC3N clumps represent long-lived, stable clumps: gravitationally bound against 
dissipation, but not collapsing. If so, then some of these regions should also be associated with enhanced 
N2H + abundances, as we in fact see in L1512. 
We now briefly consider each source in turn: 

L1512 

The strong CH3OH and SO clump in the northwest must have an age of ~ 10 4 yr. In contrast, the 
SO cores toward the southeast have only weak methanol emission and must be older. This SO emission 
surrounds the HC3N peak which must be even older, as no SO is seen at this position, suggesting an 
age of ~ 2 x 10 5 yr. This is consistent with the CO depletion and N2H+ emission also seen toward this 
position [57] 

TMC-1C 

The methanol clumps in this source have only weak SO emission, and so must be t< 5000 yr. The clumpy 
ring of material to the southeast has little or no methanol, but SO and HC3N are somewhat co-mingled. 
This is best explained by type I clumps, implying that the mechanism which formed the clumpy ring was 
not energetic enough to cause significant desorption of grain mantles. 

L1262 

In this source, the channel velocity maps at 4.12 km s _1 show a clumpy ring of SO emission surrounding 
the HC3N peak. As is the case in L1512 these may represent evolved type III clumps which have lost the 
original methanol. The central HC3N clump will be older (t> 10 5 yr) than the surrounding SO clumps 
(3 x 10 4 < t < 10 5 yr). Both CO clumps in this source display the kinematical sequence of SO followed 
by SO plus CH3OH, followed by only CH3OH, that is typical of the time evolution of type II clumps. 
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Again, we see a transition in velocity space from pure SO to pure CH3OH via mixed gas. As in L1262, 
this indicates the evolution of clumps where CH3OH ice is liberated along the line of sight. The HC3N 
peak toward the northeast is associated with weak SO, and may represent evolved gas. However, the 
peak does not appear to be correlated with the N 2 H + peak or with significant CO depletion. 

L1389 

Again, we see SO-rich gas surrounding a HCsN-rich region. As in L1512, this suggests an older, stable 
clump in the center surrounded by regions in which clump formation and/or ice mantle sublimation 
occurred more recently. This source also shows a clear HC3N-CH3OH anti-correlation. 

L1251E 

SO is only visible in this map toward the north and east, and appears to be associated with some CH3OH, 
implying an age for this gas of ~ 10 4 yr. Towards the south and west, CH3OH and HC3N clumps are 
prevalent yet SO is absent, which requires the presence of both young and old cores. This is the only 
area in which we see co-mingling of HC3N and CH3OH emission - although the peak positions are offset 
- and is the only region whose chemistry is not readily explained by our simple model. 

5 Conclusions 

We have mapped several prestellar and protostellar cores in CH3OH, HC3N, SO and C 18 0. We find that 
the emission from these molecules traces clumpy structure across all of the cores. The molecular clumps 
show striking differences in the location and velocity of emission peaks between molecules. 

We see a similar degree of chemical differentiation in all of the cores that we have observed, despite 
the differences in known star forming activity. The HC3N clumps are generally anti-correlated with 
N2H+ peaks, and with CH3OH clumps. In L1512, a HC3N clump does appear along with a N2H+ 
peak, perhaps indicative of ongoing molecular depletion. CH 3 OH clumps and SO clumps also show 
distinct emission peaks. C 18 clump peaks are not generally correlated with emission peaks in any of the 
other molecules that we observed. These morphological differences in molecular peaks towards sources 
of different evolutionary stage suggest that depletion is important, but that other processes must also be 
driving the chemical differentiation that we see. 

Our observations suggest there is a kinematic sequence through the clumps of different molecular 
species within the core. There is a well defined emission sequence in velocity space where we see only 
clumps of SO, then clumps where SO and CH3OH are intermingled, and then a region with CH3OH only. 
This may indicate that we are looking down a cylinder, or filament. 

We have investigated the origin of the chemical differentiation using a simple chemical model, in 
which clumps form rapidly and in which sublimation of ice mantles can occur as the clump is formed. 
Despite its simplicity, this model is able to account for a wide variety of the chemical differentiation 
observed in our sources. For clumps in which CH3OH is injected, the methanol survives for ~ 10 4 yr, 
so the observed CH^On-rich clumps are young compared to the time-scale for dynamical evolution. 
In these clumps, HC3N is not produced until all the CH3OH has been destroyed, which explains the 
observed anti-correlation of methanol and cyanoacctylcnc. The model also predicts some some regions 
will have HC3N/SO > 1 whereas others will have HC3N/SO < 1, depending on the age of the source 
and/or the degree of mantle sublimation. HC3N peaks at late times, and the observed regions with only 
HC3N emission are most likely tracing older clumps which are dynamically stable against both collapse 
and dissipation. Such regions should also be associated with marked CO depletion and elevated N2H + 
abundances. 

In conclusion, it appears that molecular clouds containing protostars at different epochs of star for- 
mation exhibit similar trends in their chemical spatial differentiation. The origin of this differentiation 
can be understood through simple models of the molecular gas-grain interaction, and suggests that peri- 
odic removal of ice mantles is occurring in dark clouds. Mantle disruption and dark cloud chemistry are 
probably closely connected to the physics of clump evolution, dynamical-chemical models incorporating 
of these processes are currently being developed [59] ■ 
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